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ABSTRACT 


The paper deals with the numerical simulation of the flow of Newtonian fluid by 
the control volume method, which simulates the flow of blood in a blood vessel 
with stenosis. An interpenetrating model of two-phase media is used to describe 
the process. Based on this model, under appropriate simplifying assumptions, we 
obtain the Navier-Stokes equations in the free (outside the stenosis region) zone, 
and in the stenotic region, we obtain a system of equations generalizing the 
filtration equation. Stenosis is viewed as a porous medium. The influence of the 
Reynolds number and the degree of stenosis on the nature of the flow on the 


behavior of blood is investigated. 
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INTRODUCTION 


Arterial stenosis occurs due to the deposition of cholesterol on its walls. Stenoses 
restrict blood flow. This change in blood flow caused by stenosis can lead to the 
development of arterial disease and arterial deformity. 

With stenosis of the artery, the nature of the blood flow in the lower part of 
the artery changes significantly. These changes are highly dependent on the 
degree of stenosis. Due to the presence of moderate or severe stenosis in the lower 
course, the course is severely disrupted and, therefore, the course becomes 
unevenly complex. 

Many studies have been investigated to understand the nature of blood flow 
disturbances in stenosis. Both experimental and numerical methods have been 
used by various researchers to determine the development of stenosis and its 
magnitude [1-4]. 

It is assumed here that the movement in the vessels is laminar. To simulate 
the flow of blood in vessels with stenosis, we use an interpenetrating model. We 
consider stenosis as the second phase in an interpenetrating model for describing 
heterogeneous media with a fixed zero velocity. The area free from stenosis, the 
blood flow is described by the Navier-Stokes equations. We represent stenosis as 


a porous medium. 


Flow description model 
Consider an interpenetrating model describing the flows of two-phase media in a 
Cartesian coordinate system (two-dimensional problem) [6-7]: 
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Here, us — the j-th velocity component, p — pressure, 7 — volumetric concentration of the 1 — phase ee = 2, / Poi), ro) in 
is the Kronecker symbol, {/; — the viscosity of the 1— phase, K —the phase interaction coefficient, O,, Qo; — are the reduced and 


true density of the i— phase ( Pp; — const), p? — the j-th component of the 7 — phase mass force. 


If we represent in the equations that the second phase is stationary, we obtain a system of equations describing the flow of a liquid 


through a stationary porous layer: 
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As accepted in porous media, we take the interaction coefficient in the form: 


_ 150(1- f)’ 
= ar 


Here d — is the characteristic size of the porous layer. 


It is easy to see that at f =1 equations (4), (5) transforms into the Navier-Stokes equations. When 14 f <1 we get a system 


of filtration equations. Let us write equations (4) - (5) in dimensionless form (stationary case): 
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where u, v— are the longitudinal and transverse components of the velocity, x, y — are the Cartesian coordinates, Re = hUp/u— is 


the Reynolds number, U— is the characteristic velocity, h — is the characteristic size, D=~vahild is the dimensionless number, 
d — is the characteristic size of the porous medium, a — coefficient of proportionality. 


For the numerical solution of the system of equations (6) - (8), the control volume method [8] with appropriate generalizations was 
applied. An uneven coordinated mesh is built by thickening near the stenosis. 


BOUNDARY CONDITIONS 


At the entrance to the vessel, a parabolic law of the velocity distribution is given, at outputs a soft condition, on the walls of the 
vessel, the no-slip condition: 


x=0, u=6y(l-y), v=0, p=p’. 


With sufficient distance, we set the soft condition 


Y=. a au 


== 0. 
Ox Ox 


STENOSIS MODELING 


Consider to determine the contribution of stenosis in blood vessels based on the model (6-8). As a form of stenosis, let us take the 
form given in [9]. We believe that the stenosis is asymmetrical (Fig. 1). The form of stenosis is determined by the equation 


y(x) = Mo Pecos er ee 
es L, ps 


Page2 2 / 


ISSN 2278-5469 EISSN 2278-5450 | OPEN ACCESS 


DISCOVERY | REPORT 


Figure: 1. Area of flow with stenosis 


Here d is the distance from the entrance to the vessel, hy is the fullness of the vessel, L) is the width of the stenosis 


NUMERICAL EXPERIMENTS 


In the calculations, the parameters of stenosis were taken: d =0,5; Ll, =1. The parameter hy varied. The calculation was carried 


out in the area[O,1]x[0,L] (L = 4). Number of nodes 50x30. 


05 0,75 1 1,25 15 1,75 2 2,25 25 2,75 3 3,25 3,5 3,75 


Figure: 2 Velocity distribution over the sections of the vessel, the fullness of the vessel is 20% 


In fig. 2 shows the distribution of the velocity in different sections: the asterisks represent the distribution of the velocity in the 
section: diamond - in cross section x = 0,45 ; rectangles - x = 1,025 and circles - x = 3,337. The crosses represent the border 


of the stenosis. 


05 0,75 1 1,25 1,5 1,75 2 2,25 25 2,75 3 3,25 3,5 3,75 


Figure: 3 Velocity distribution over the sections of the vessel, the fullness of the vessel is 40% 


Figure 3 shows the calculation results for the same parameter values. The difference lies in the fullness of the vessel. In this case, 


there is a circulating flow behind the stenosis in the area X € [1,3;2,56]. For example, in Fig. 3 the graph of the velocity 
distribution corresponding to the cross section x =1,809 in the lower part is observed reverse flow. However, with the same 


parameters, but with porosity f = 0,6, there is no circulation flow. 
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With an increase in the filling of the vessel (Fig. 4), an increase in the zone of reverse flows is observed. The secondary flow 


covers the zone x €[1, 2533, 8]. 


Figure: 4 Velocity distribution over the cross-sections of the vessel, the fullness of the vessel is 60% 


CONCLUSION 


On the basis of an interpenetrating model of the flow of heterogeneous media, a mathematical description of the flow of blood with 


porous stenosis is proposed. For two-dimensional problems, a discrete model of the equation is constructed using the control 


volume method. The results of calculations and the influence of grid parameters, reflecting the legitimacy of the application of the 


model and the proposed numerical approach, are presented. 
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